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Abstract 



In this paper, we present the detailed calculation of the persistence exponent 
8 for a nearly-Markovian Gaussian process X(t), a problem initially intro- 
duced in [Phys. Rev. Lett. 77, 1420 (1996)], describing the probability 
that the walker never crosses the origin. Resummed perturbative and non- 
perturbative expressions for 6 are derived, which suggest a connection with 
the result of the alternative independent interval approximation (HA). The 
perturbation theory is extended to the calculation of for non-Gaussian pro- 
cesses, by making a strong connection between the problem of persistence 
and the calculation of the energy eigenfunctions of a quantum mechanical 
problem. Finally, we give perturbative and non-perturbative expressions for 
the persistence exponent 6(Xq), describing the probability that the process 
remains bigger than Xq x yj (X 2 (t)). 



PACS numbers: 02.50, 05.40.+j, 05.50.+q, 82.20.Fd 
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I. INTRODUCTION 



A natural quantity that characterizes a given stochastic process X(t) is its persistence 
P(t), i.e., the probability that this signal has kept the same sign up to time t. For a 
large class of physical systems (to be defined more precisely below), persistence decays as a 
power-law in time, P(t) ~ t~ 8 for large t, thus defining the persistence exponent 9. 

This exponent has been studied in experimental systems (breath figures Q, a liquid 
crystal system mimicking the 2d Ising model 0, soap bubbles 0]...), an d by theoretical 
means through the exact solution of models [§-0], numerical simulations f||9[], and general 
theoretical methods JiTf-fTo]] . 

Most theoretical methods restrict themselves to the study of persistence of stochastic 
processes that are Gaussian. This is partly because Gaussian processes are abundant and 
simpler. Moreover, in many physical situations, the study of persistence of non-Gaussian 
signals can be effectively reduced to that of Gaussian signals [0,0 • Thus, given a Gaussian 
process X(t) of zero mean, the basic question is: what is the probability that it remains, say, 
positive up to time t? This is a difficult problem that has been studied by mathematicians 
for a long time |I6| |. Recently, however, it has created much interest among physicists. 

One of the general methods recently introduced to tackle this difficult problem, namely 
the independent interval approximation (IIA) ||12|| , assumes that the interval lengths be- 
tween zeros of the process X(t) are statistically independent. This sole assumption permits 
the closure of a hierarchy of equations leading to an approximate expression of 9. This 
approximation gives very good results for smooth Gaussian processes {i.e. processes with a 
finite density of zero crossings). Unfortunately, it is not clear how this assumption can be 
justified, and whether a Gaussian process can be said a priori to be well described by this 
approximation. 



An approximation for the distribution of the time-integrated "magnetization" fli"3| . |T4 
M(t) = t" 1 / sign(X(u)) du, also leads to good quantitative results for smooth processes, 
but suffers from the same conceptual problems as the IIA, and is not even guided by physical 
intuition. Still, the study of this quantity has lead to the introduction of a new quantity, 
the generalized persistence [|13|| , that is the probability that M(t) remains above a certain 
level M . This quantity which decays with a persistence exponent depending continuously 
on Mo has been studied in the framework of spin systems and for random walkers JIB| . 

Finally, a systematic e-expansion, which is exact order by order, has been developed 
recently for smooth Gaussian processes 

However, all these approximate and exact techniques fail for processes that are singular, 
that is for which the density of zero crossings is infinite. These processes appear in many 
physical situations such as nearly Markovian random walkers or interface growth [O. 

In this paper, we come back to the first general method proposed, that is perturbation 
theory around a Gaussian and Markovian process [TD|. After introducing the principle 
of this method (section II to IV), which shows a deep connection between the problems 
of persistence and the energy spectrum of a quantum mechanical problem, we present a 
symmetry argument for the exact form of 9 which leads to more general results for the 
persistence exponent (section V). These results also reveal a connection between IIA and 
perturbative approaches. In section VI, we extend the perturbative approach to the case 
of non-Gaussian processes further reinforcing the link with standard quantum mechanics. 
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In section VII, we show that the various approaches introduced can be applied to the 
computation of the probability that the signal X(t) / yj (X 2 (t)) remains higher than any 
given non-zero-level Xq (generalized persistence). Finally, in section VIII, we illustrate 
some of the results obtained in the preceding sections by means of numerical simulations. 



II. IMPORTANCE OF GAUSSIAN STATIONARY PROCESSES 



The most popular examples of persistent systems have been taken from the field of 
coarsening dynamics [[L8[]. For instance, let us consider an Ising spin system after a quench 
at very low temperature from a high temperature disordered state. Domains of positive 
(essentially +1) and negative (essentially —1) magnetization grow with a time-dependent 
typical length scale L(t) ~ t 1 ! 2 . For this system, the spin persistence, that is, the probability 
that a spin has never changed sign, or has never been crossed by an interface, is known to 
decay as t~\ with 9 = 3/8 in d = 1 @, and 9 « 0.22 in d = 2 @-|nJ. 

Due to dynamical scaling the two-time spin correlation function only depends on the 
dimensionless ratio of L at both considered times: 



(S(t)S(t')) = f(L(t)/L(t>)). 



(1) 



This property will be characteristic of a coarsening system and only relies on the existence 
of a unique dynamical length scale and the dynamical scaling hypothesis [JTB|. Now, if 



L(t) behaves as a power-law for large times, all two-point correlation functions are then 
functions of t/t'. By considering r = log(t), these correlation functions are then functions 
of exp \t — t'\, or more simply |r — r'\, so that they become stationary in the fictitious time 

T. 

Moreover, in many physical systems |T0| , p!2|JT9[| , the question of computing the persistence 
for the original dynamical variable (S(t) in the above example) can be reduced to the 
study of the persistence of a Gaussian variable X(t). One possibility is, of course, that the 
physical variable is a Gaussian variable itself: this occurs in the study of the persistence 
of the diffusion equation [12], and for the total magnetization persistence of a spin system 
quenched at T < T c [p0[] , or T = T c (in the latter case, the persistence exponent is a new 
critical exponent |19|). But in some other cases, including the Ising and more generally 



O(N) spin systems, the original persistence can be shown to be very close to that of a true 
Gaussian process X(t). For instance, a local spin in an Ising system behaves essentially 
as the sign of such a Gaussian process, S(t) « sign(X(t)), an important result which was 
first used within the OJK theory pT,18|, and later more precisely formalized by Mazenko 
and co-workers 



22,18 



To summarize, we underline the special role played by Gaussian 
processes, and will thus restrict our study to this kind of process. 

The next important remark is that if the persistence of the considered Gaussian process 
X(t) decays as t~ 9 , the persistence in terms of the fictitious time r (for which this process 
is stationary) is expected to decay exponentially as exp(— Or). Thus, in the following we 
restrict ourselves to the study of persistence for a stationary Gaussian process X{r) p3|-|26|. 
Note that if L(t) does not behave as a power-law of time, the persistence still decays as a 
power-law of L(t) as soon as dynamical scaling is satisfied, and the proper fictitious time is 
simply r = logL(t), for which the process X(r) is again stationary. 
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The most general equation of motion for a stationary Gaussian walker reads, 



X'(r) = -XX(r) + f J(r - r'Hr')dr' , (2) 

J — oo 

where t](t) is a Gaussian white noise satisfying (v(' r )v( T> )) = S(t — t'). Indeed, this equation 
must be linear to preserve the Gaussian property, and the coefficient A of X(r) must be 
constant to preserve stationarity. The last term of Eq. (|2]) accounts for memory effects, 
involving a memory kernel J, and must take the form of a convolution product, again to 
preserve stationarity and the Gaussian property (linearity). Note that it is not necessary to 
involve higher derivatives of X in this equation of motion, as they can be accounted for by 
a proper choice of the kernel J (see Eq. (EJ) below). 

The Markovian case is associated with J(r) = 5(r) (no memory effects), so that the 
equation of motion becomes, 

X'(r) = -XX(r)+r ] (r). (3) 

The velocity X'(r) only involves the noise at the same time r. For such a Langevin walker, 
the two-point correlation function is simply, 

(X(r)X(r')) = f(T-T'), /(r) = eXP 2 ~ A|T| . (4) 

For convenience, the correlator (and the variable X) has been normalized such that f'(0 ± ) = 
=pl/2, and, from now on, this will be assumed for all correlators. This will ensure that, 

oj 2 f(uu) — > 1, when uo — > ±oo. (5) 

Also note that this Markovian correlator / has a cusp at the origin. We will define a nearly 
Markovian Gaussian process as one with a correlator which satisfy the above condition (5). 

In general, the knowledge of the two-point correlation function /(r) is equivalent to that 
of the equation of motion, as the Fourier transform of / satisfies, 



f{u) = (X(co)X(-co)) = J^t. (6) 



This actually shows that any correlator f(ui) can be reproduced by a proper (not unique) 
choice of the memory kernel J. 

In sections III and IV, we will give a more extensive account of the perturbative ex- 
pansion for 9, in the case of a nearly Markovian Gaussian stationary process, a calculation 



which was first introduced in ||1Q|| , and then reproduced in a real time formalism in [ IT 
This will be followed (section IV) by a resummation of this perturbation theory using a 
general symmetry argument, and the discovery of an intimate connection between the IIA 
and perturbative methods. A new non-perturbative expression for 6 is also presented, which 
happens to reproduce quantitatively most numerical results (section VIII). 
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III. PERSISTENCE: THE MARKOVIAN CASE 



Let us now move to the problem of persistence. The probability that a given walker 
remains on, say, the positive side of at all times between and (3 is, 

f x>0 VX(r)eM-S] 2h 
nP) JVX(r)exp[-S] Z > { ' ) 

where, 

509, {X(t)}) = i /* IfnMn-r^jdrA, (8) 

is the Gaussian weight associated with the trajectory X(r), and g(j\ — T2) is the inverse of 
the correlation matrix /(n — T2). 9 is then calculated from P(/3) by taking the limit, 

B = - lim /TMogP^). (9) 

We can impose periodic boundary conditions for the walker trajectories, X(0) = X({3), 
which should not affect the value of 6 in the limit of large (3. Indeed, in practice, the process 
will have a finite typical correlation time, equal to A -1 in the example of the Markovian 
walker, so that this extra constraint cannot affect the large time persistence regime. 

The path integrals of Eq. ([?]) strongly suggest the connection of this problem to Feynmann 
integrals in quantum mechanics or statistical field theory. Let us make this connection more 
precise. Because of the periodicity of the trajectories, the Gaussian weight in Eq. (|7|) can be 
also written, 

1 +00 

5 = ^EsK)|x(^)l 2 , (10) 

Z P n=0 

where g{oJ n ) = l/f(u n ) (the kernel in the expression of S is diagonal in Fourier space) 
and u n = 2irn/(3 are Matsubara frequencies. First consider a Markovian process for which 
g{ui) = uj 2 + A 2 (the Fourier transform of /(r) = exp(— A|r|)/2A is [uj 2 + A 2 ] -1 ). S can be 
alternatively written as, 




dr. (11) 



We recognize the action in imaginary time (f3 is then an inverse temperature) of an harmonic 
oscillator of frequency A. The periodicity of the paths ensures that Z =Tr [exp — /3H ] 
is then the partition function of an harmonic oscillator, and Z\ =Tr [exp — /3i?i], is the 
partition function of the same harmonic oscillator, but with an infinite wall at the origin (as 
the particle is constrained to remain on the positive axis). For large time, the persistence 
behaves as, 

P^-exph^-i^)], (12) 

where E and Ei are the ground state energies of these quantum systems. By direct identi- 
fication, we thus find that, 



5 



(13) 



Moreover, En = A/2, and it is easy to convince oneself that the ground state wavefunction 
of Hi, is the first excited state of H restricted to the positive axis, so that Ei = 3 A/2 (this 
argument is very general, and only relies on the x — * —x symmetry of the potential). We 
finally find that 9 = A for a Markovian process. This is a well-known fact [^ -|26|, that can 
be simply illustrated for the usual Langevin Markovian walker, for which the equation of 
motion reads (in actual time t), 4| = f]{t)- For such a random walk, the persistence exponent 
is known to be 1/2 p3|-p6|. Let us reproduce this result within our approach. The two-point 
correlation function is easily computed: (x(t)x(t')) = min(i, £'), and the normalized variable 



X(t) = x{t)/J{x{t) 2 ), has a correlator, (X(t)X(t')} = (t'/t) 1/2 , for t > t'. This correlator is 
a function of the ratio of the two times, so that it is stationary after the change of variable 
t = log(t), becoming (X(t)X(t')) = exp —\\t — r'| . Applying the above calculation, we 
recover the result, 6 = A = 1/2. 



IV. PERTURBATION AROUND A GAUSSIAN MARKOVIAN PROCESS 

Of course, this heavy machinery is not introduced to deal with the well understood 
Markovian case, but rather to be applied to the case of a nearly Markovian walker, for 
which no result exists. Thus, let us consider such a walker for which, 

f(r) = ^ [exp(-AM) + 0(r)] , (14) 

where 0(r) is assumed to be a "small perturbation" to the Markovian correlator. In Fourier 
space this can be written, to first order in 0, 

g{u) = f{uo)- 1 = uj 2 + A 2 - h{u), h(cj) = ^g/^ M- (15) 

In the general case, the denominator Zq of Eq. (|7|) can be exactly computed, as any uncon- 
strained Gaussian integral, and is proportional to det 1 ^ 2 [/(r i —Tj)\. After taking the proper 
limit, Eq = — lim^+oo (3' 1 logZ (/5), one finds, 

E = -— log L 2 f (cu)) dco. (16) 
Zn Jo v 7 

Note that this integral converges thanks to the relation expressed in Eq. (H). To be consistent 
with the perturbative expansion for E\ to come, we can write i?o up to first order in <f>, 

E ° =2-^xJ (^ 2 + A 2 )0M^ + O(0 2 ), (17) 

the first term being the previously discussed Markovian result, that is the ground state 
energy of an harmonic oscillator of frequency A. The computation of Zi (or E{) is still a 
formidable task, as the domain of integration of the Gaussian integral only involves positive 
X(t), for all r. The natural impulse is to write S = S osc , + 5S, where S osc . is the harmonic 
oscillator action associated with a Markovian process (Eq. (|HD), and, 
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SS = - I fg jg XinWn - r 2 )X(r 2 ) dndn, (18) 
= -^£^AK)|*K)| 2 , (19) 

where the Fourier transform of h is given in Eq. (|1|) . 5S is linear in <p and can be considered 
as a small perturbation to <S osc .. We can now use the standard first order cumulant expansion 
of quantum mechanics (or statistical field theory), leading to, 

E l = — + lim (8S) wall + O(0 2 ), (20) 

2 /3-»+oo 

where the average is to be taken using the Boltzmann weight associated with the harmonic 
oscillator of frequency A, with an infinite wall at the origin. Let us denote by \T) the 
eigenstates of this quantum system (as opposed to |Z), the eigenstates of the unconstrained 
oscillator), associated with the eigenenergies ei = {{21 + 1) + 1/2)A = {21 + 3/2)A (/ > 0). 
One can then write, 

(0\X(-u n )X(u n )\0) = / 2cos^r^|(0|X|/)| 2 e- (£i - £o) ^r. (21) 

Jo 1=0 

|(6|X|/)| 2 can be computed for the harmonic oscillator with a wall, using the fact that 
(x\l) = y/2(x\2l + 1), for x > 0, and exploiting standard properties of Hermite polynomials. 
The complete calculation is performed in appendix A and B. The final result reads, 

( 6|x(-.„)x(.^6> = A,(^) + g_||_, (22) 

the Dirac peak coming from the I = term. The coefficients Cj involved in this relation 

read, 

- i ( V ( 23 ) 

Finally, the sum over n in 8S becomes an integral in the (3 — > +oo limit, leading to the final 
expression for 9 = E\ — Eq: 

I r+oo A 

9 = X-— V{u)(f){u;)du + 0{(f) 2 ). (24) 

27T JO 

The kernel V is defined by, 



-6 (w) + ^ 



A ^w 2 + 4j 2 A 2 u 2 + X 2 



(25) 



As noticed by Oerding et ai, this cumbersome expression in terms of Fourier transforms 
has a remarkably compact form when expressed in the inverse Fourier space. Indeed, the 
function between brackets is just the Fourier transform of, 

I +oo 1 

U{r) = -J2c J exp{-2jX\r\) - — exp(-A|r|), (26) 
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with c = -, so that V is the Fourier transform of 4\(—d 2 + \ 2 ) 2 U{r). This allows us to 
recast the preceding result into the form, 



e = a 



2A Jo 



r)(-d 2 T + X 2 yU(r) dr + O^ 2 ). 



(27) 



A simple manipulation on the c/s (see appendix A) allows to resum exactly the series 
(-d 2 + \ 2 ) 2 U(t), finally leading to, 



e = a 



1 _ _ / 0( r )[i _ exp(-2Ar)]- 3 / 2 rfr 

7T JO 



O(0 2 



(28) 



We can generalize this expression when the constraint on X(r) is X(r) > X , instead 
of X(t) > [[T3|j27[1 . Indeed, for the Brownian walker (f(r) = exp(— |r|/2), such that 
(X 2 ) = /(0) = 1), it is known (see section VII) that 9 satisfies V 2e (X ) =0^, where 
D28 is a parabolic cylinder function. We can expand this expression for small Xq, leading 
t° ^Brownian = 1/2 + XqJ \^2n + O(X ) 2 . If we perturb around a general Markovian process 
(/( r ) = ix ex P( — ^l r l))> we then get another perturbative contribution for the exponent 9 
(valid in the limit of small Xq), which should be added to the result of Eq. (|28|): 



89(X 



2A x 



2tt(X 2 ) 



ow 2 



f = A 3/2 • 



2X 



7T 



O(X ) 2 . 



(29) 



V. RESUMMATION: A SYMMETRY ARGUMENT 

A. Resummation in time 

Consider a Gaussian process of correlator / and persistence exponent 9. Let us assume 
that we have been able to resum all terms of the perturbative expansion which contain only 
one time integral. Very generally, one can thus write, 

r+oo 

6= A(f(r)/f(0),r)dr. (30) 
Jo 

The variable /(t)//(0) appears due to the fact that 9 should not depend on the correlator 
normalization (here /'(0 ) = but /(0) = 1 was chosen in 0,0). 

If /(t) is changed into /(err), it is clear that the persistence exponent is simply changed 
into ot6. Using this remark, we get, 

r+oo 

a6= A(f(ar)/f(0),r)dr, (31) 

J 

which shows after a simple change of variable that, for any process and any a, one must 
have, 

r+oo 

9= A(f(r)/f(0),r/a)dr/a 2 . (32) 
Jo 
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This strongly suggests that 6 can in fact be written as, 



f+oo A T 

0=1 B(f(r)/f(0))-. (33) 



r 2 



Assuming now that /(r) is close to a Markovian process with an associated small 0(t), one 
can develop Eq. (^) leading to, 

6 = / £?(exp(-Ar)) % + / 0(r)S'(exp(-Ar)) - + O(0 2 ). (34) 
jo r A Jo T A 

In this perturbative limit, Eq. ( |34"D should coincide with Eq. (^) leading to, 



B'(exp -X) = X \ (35) 

7r (1 — exp — za)^ 

or, after making the change of variable u = exp —X, 

2 \og(u) 2 

B M = -n(l-u*)W (36) 

For the integral Eq. (|33| ) to converge one should have B(l) = 0, finally leading to the final 
result of Eq. ([33]), with B given by: 

*(«> = !/' 7=^*. (37) 



7T 



« (1-W 2 ) 3 / 2 



Note that this expression is not only defined for a nearly Markovian process, for which / 
has a cusp at r = 0, but actually converges for any process for which, 

/(r)//(0) - 1 ~ |t|", when r - 0, (38) 

for any /x > 2/3 (as 5(1 — e) ~ e 3 ^ 2 ). Smooth processes (with a continuous velocity) are 
associated with fi = 2, and the local density of a charge distribution evolving according to 
the simple diffusion (or heat) equation corresponds to this case [Q. As we have argued 
in detail, fi — 1 corresponds to nearly Markovian processes. Finally, other values of \i < 2 
correspond to singular walkers for which the fractal density of the set of X = crossing 
times is 1 — /i/2. Such processes have been encountered in the study of out of equilibrium 
atomic surfaces, for which X(t) is the local height of the substrate [ |i~7| . 

A nice consistency check consists in showing that the first term in Eq. fl3~4] ) is equal to 
A, that is the Markovian value for 9. This is simply done by performing an integration by 
parts using the explicit expression of B, leading to, 

2A /•+ oc (Ar) 2 e- Ar dr 2A f 1 logw 2A f 1 du 

C'Markov — / TTo ' — / oTo du — / — , — A, (39) 

TV Jo (1 - e- 2Ar ) 3/2 T 7T Jo (1 - u 2 ) 3/2 TT Jo y/1 — U V ' 

the last integral being obtained through another integration by parts. 

The argument presented above was motivated by the following important remark: for a 
given correlator /, the perturbation <fi, or equivalently, the function exp(— A|r|)/2A around 
which the perturbation is started, are actually quite ill-defined. If we knew the complete 
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perturbation expansion, starting from any value of A we should get the same result. Note 
that in standard field theory, one usually perturbs around a system which is solvable for a 
certain value (usually 0) of the coupling constant: there is a unique way of performing the 
perturbative expansion. Thus, it is natural to ask whether there is an optimal choice for 
the starting value of A. A very natural choice is to take for A the value which cancels the 
first order perturbative term. In other words, we take the "best" starting Markovian cor- 
relator such that the first order contribution vanishes. This gives another non-perturbative 
expression for 9 (that we may call "variational" or self-consistent perturbative), which must 
satisfy, 

/W//(0)-exp(- e ,) (ir = a 



[1 - exp(-2#r)] 3 / 2 

This equation always has a solution, as the expression in Eq. ( |£0| ) is clearly positive for 
9 — > +oo, and goes to — oo, when 9 — > 0. The expression in Eq. (|40|) is defined for any 
\i > 1/2, in fact a larger domain than the fully resummed formula of Eq. fl3"7p. 
Note that we can write the resummed expression in a similar form, 

/ [B(/(r)//(0)) - B(exp(-9r)} dr = 0, (41) 
Jo 

which after integration by parts, takes the form, 

-~tf[/(r)//(O),exp(0r)] 



o (l-/ 2 (r)// 2 (0)) 3/2 



dr = 0, (42) 



where the precise form of the known kernel K is of no real interest. This last remark allows 
us to make a link with the IIA result. Within this scheme, based on the approximation that 
the intervals between the zeros of the process are independent, it can be shown for smooth 
processes (// = 2) that 9 must satisfy Jl2| (with the normalization /(0) = 1), 



tx9 



29 r+°° 

1 + — / exp(#r) sin^f/fr)) dr 
n Jo 



-no) 

If one integrates by parts this expression twice, it takes the following form: 



0. (43) 



+oo 







This expression now looks of the same type as the ones found within the perturbative 
approach. However, its domain of definition remains strictly — 

Finally, let us mention that when the constraint on X(r) is X(r) > X [|13|,^7]] instead of 
X{t) > 0, the following perturbative correction should be added to the preceding expressions 
for 9 (see Eq. 

69(X ) = 9(X = 0) • ==L= + O(X ) 2 . (45) 

/2tt/(0) 
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B. Resummation in frequency space 

The same argument as above can be applied to the expression of 9 in frequency space. 
This time, this will allow us to resum all terms in the perturbation theory involving only 
one frequency integration. Again, we assume that 9 can be written, 

r+oo 

9=1 C{f{u),u)du. (46) 

J 

We still assume that / has a finite derivative in + , keeping the normalization 2|/'(0 + )| = 1. 

If /(r) is changed into f(ar)/a (to preserve the normalization), f(u) is changed into 
f(co/a)/a 2 , and it is again clear that the persistence exponent is simply changed into a>9. 
Using this remark, we get: 

aB = J C [a- 2 f(u/a),u) du, (47) 

which shows after a simple change of variable that, for any process and any a, one must 
have, 

r+oo . . 
9= C (a- 2 f(u),auo) duo. (48) 

J 

This again strongly suggests that 9 can in fact be written as, 

9= / D{uo 2 f{uo))duo. (49) 
Jo 

Note that this property is shared by the exact and general expression of Eq , given in Eq. ([U]) . 

Now, one can consider a nearly Markovian process, for which the correlator satisfies 
Eq. (|l"5|). One can develop Eq. (JEJ) up to first order in <f) and identify the result to the 
perturbation result of Eq. fl24j - |25D . The calculation is elementary and leads to, 

4 - 1 r+oo A 

6 = -^(O) + — / W(u 2 f(uj)) duo, (50) 
7r Ztt Jo 

+oo _ 

W(x) = ^ — log (l + An 2 (x- 1 - l)) + log(x). (51) 

n=l 71 

The first term arises from the 5 term in the kernel V, and can also be written as 
f Jo +0 ° ^ ( y UJ \f( UJ )\ 1 ^ S ) du. Again, it is easy using relations given in appendix A to check that 
the Markovian value 9 = X is recovered for the Markovian correlator f{uo) = (u 2 + A 2 ) -1 . 
Note finally that this procedure permits the recovery of the exact expression of Eq, which 
is produced by the last log(x) term in the kernel W. 

VI. PERTURBATION AROUND A NON-GAUSSIAN MARKOVIAN PROCESS 

When writing S = S osc , + SS, we deliberately chose to perturb around a Gaussian 
Markovian walker, or around the quantum action of an harmonic oscillator in terms of the 
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pseudo-action S. This is quite arbitrary, and in principle the action of any (preferably 
solvable) quantum system would have worked. The stochastic process associated to such 
an action (each trajectory {X(r)} being weighted by exp — «S[{X(r)}]) is Markovian but, 
in general, non-Gaussian, as the only Gaussian quantum action is that of an harmonic 
oscillator. 

So in this section, we consider a stationary stochastic process X(t) of any kind, associated 
with the weight or pseudo-action S, and a quantum mechanical system for which the action 
is Sq ("Q" for "quantum"). This quantum system could be an harmonic oscillator, a particle 
in a square box, and more widely, any system preferably solvable for the actual perturbative 
calculation to be tractable. 

Then, setting S = Sq + SS, and reproducing exactly the calculations of the beginning of 
section IV, we end up with the following perturbative expression for 9: 

9 = E?-E$ + lim [{SS^-iSSU + OiSS 2 ), (52) 

f3— >+oo 

where (respectively E®) is the ground state energy of the unconstrained (respectively 
with an infinite wall at X = 0) quantum system. ( )o and ( )i denote quantum averages 
performed using the Hamiltonian of the quantum system, respectively without and with the 
infinite barrier at the origin. 

We have already implicitly used Eq. fl52|) in section IV, where Sq was chosen to be the 
Gaussian quantum action of an harmonic oscillator. Let us now illustrate Eq. ([52]) by taking a 
non-Gaussian system as the starting quantum system. The simplest possible example is that 
of particle in a box, with X restrained to the interval [—b, b] . We now use this simple non- 
Gaussian system to compute approximately the value of 9 for a Gaussian process associated 
with the Gaussian weight S defined in Eq. (||). 

Let us call \l) (I > 0), the eigenstates of a quantum particle in the box [—6, b], associated 
with the eigenenergies E\ — \k 2 (l + l) 2 , with k = ^. The eigenstates of the constrained 
system (a particle in the box [0, b}) are the \l) = \/2|2Z+l) (/ > 0). To evaluate (5S)i — (SS) , 
one essentially needs to compute (0|X(— u n )X(uj n )\0) and (0\X (— u> n )X(uj n )\0) . This is a 
straightforward task using identities similar to Eq. (|2l|), where the scalar products, (0|X|/) 
and (0|X|/) are even easier to compute for a particle in a box (see appendix B). Introducing, 

95fi +°° n .ft. 

fcK) = (0\X(-u n )X(u n )\0) = — £ ' J , , (53) 



with, 



7t 2 Kb] + ^1 



2f?' + Ll 2 1 

and bj = -(2j + 1)(2j + 3), (54) 



J (2j + 1)4(2j+3) 4 ' J 2 

we get the following expressions for E 1 and E {9 = E\ — E ): 



kl 1 f +oc ( k(u) \ , 



and, 
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^ = 2fcg + ^ ( 7T 4^ - ] (56) 




32tt Jo 



fc(u;/4) dw. (57) 



For a given correlator f(u), it is not clear what the "best" starting value for k (or for 
the box size b) is. Let us propose two natural choices. We can take ko such that the first 
order perturbation vanishes, which leads to 9 = 3A;q/2. An alternative choice is to take ko 
such that Ei is minimum, as it can be shown that Ei(ko) has always such a minimum for a 
finite ko- In fact, the variational inequality E\ < 2k\ + lim^ +00 (<5iS)i is exact for any ko, 
which intuitively validates this choice if k Q . 



VII. GENERALIZED PERSISTENCE 

So far, we have essentially considered the probability that the signal X(t) has never 
changed sign. In fact, it seems natural to study the more general probability that the 
signal has always remained above a certain level X . When X ^ 0, this defines the X - 
level persistence. This generalized persistence has already been introduced for the simplest 
Markovian Gaussian walker |27| , and spin systems [HJ . Moreover, at least in the framework 



of Mazenko approximation ||22|| , there is a connection between the persistence of the q- Potts 
model RP|| (the probability that a given site always remains in a given phase) and the 



Xq = F(g)-level persistence of a certain Gaussian variable [28 |. 

Let us take the example of the Gaussian Markovian walker, associated, within our for- 
malism, with the action of an harmonic oscillator. E\{Xq) is now the ground state energy 
of an harmonic oscillator with an infinite barrier at X , for which the eigenstates can be 
expressed in terms of a parabolic cylinder function |27|] (generalization to a continuous index 
of Hermite polynomials). E\ is then implicitly defined by imposing that the ground state 
eigenfunction has a unique node at X = X . If we come back to real time t = expr (see 
section III), Xo-level persistence for the Langevin walker satisfying ^ = r)(t) is defined as 

the probability that x(t) always remained greater than Xo\J (x 2 (t)) = X Vi. This decays as 
a power-law of time with exponent 9(Xq) = Ei(Xq) — Eq. 

If we were to compute the A -level persistence exponent 9(X ) for a Gaussian process 
using the perturbation theory formalism, we would have to evaluate scalar products like 
(0|A|/), where \l) are the eigenstates of the harmonic oscillator with an infinite barrier at 
X = Xq. Unfortunately, this seems to be an analytically untractable problem. However, 
using the formalism of the preceding section, we only have to evaluate brackets involving 
eigenstates of a particle in the box [Xq, b], which are explicitly known. The calculation is 
straightforward (see appendix B) and leads to, 

Ei(Xo) = + ^ 2 (1 + V) 2 (tA^ - T^b) (58) 

k(u(l - rif/A) dw, (59) 




13 



where 77 = Xo/b, and k has been defined in Eq. (J53Q . As a check, we can see that for 77 = 
{i.e. the wall is at the origin) we recover the result of Eq. (|57|), and for 77 = —1 {i.e. the wall 
is at X = —b, which corresponds to no effective constraint), we recover the expression for 
E of Eq.(||). 

Again, for a given process X and a given level Xq, ko can be fixed by imposing that the 
first order perturbation term in 9 vanishes, or by taking the value of ko for which Ei{Xq, ko) 
is minimum. 



VIII. NUMERICAL SIMULATIONS 

We now illustrate the various analytical results obtained in the preceding sections by 
means of numerical simulations. 



A. Nearly Markovian processes 

As already mentioned in the introduction, a local Ising spin evolving after a quench 
at T = 0, from the high temperature disordered phase, essentially behaves as the sign of a 
Gaussian variable. Mazenko approximation |22| then permits the calculation of the two-time 



correlator of this Gaussian process. It happens that, in one dimension, this approximation 
recovers the exact expression of (S(t)S(t')) |18[], leading to the following form of the correlator 



/ when expressed in the fictitious time r = log(t): 



fir) = \h- r- v (60) 

V 1 + exp(r) 

The exact value of 9 in d = 1 is 6 = 3/8 = 0.375 [|J. The "variational" and resummed 
perturbative expression of Eq. ( f40~|) and Eq. (|33| - |3~7|) respectively lead to 9 v3ir = .3595 . . . 
and 6* P ert = .3677.... The process associated with the correlator given by Eq. (|60|) has 
been actually simulated using the Fourier space form of Eq. (|2]). We have obtained 9 = 
0.355 ±0.005, in extremely good agreement with the theory. The small discrepancy with the 
exact result 0.375 for the Ising model is attributed to the fact that the actual process such 
that S(t) = sign(A(t)) is not strictly Gaussian. However it seems that this non-Gaussian 
effect is rather small. 



We have also tested our theoretical expressions using a correlator introduced in [1C 



/(r) = -exp(-r) + -exp(-2r). (61) 

We found 9 vax = 1.4855 . . . and pert = 1.4802 . . ., again in good agreement with the numerical 
result 9 = 1.481 ± 0.005. 



B. Other singular processes 

Interesting examples of singular correlators with fi < 2 and /ia 7^ 1 (see the definition of 
/1 in section V) have been introduced in the framework of dynamical surfaces described by 
the following time-dependent equation |17| : 



14 



(62) 



where h is the local height of the fluctuating interface, and rj is Gaussian white noise. The 
equation being linear, h(x, t) is a Gaussian variable, for which we take the initial condition 
h(x, 0) = 0. To define the first passage problems of interest, consider the quantity, 

P(t , t) = Prob[/i(x, s) ^ h(x, t Q ) Vs : t < s < t + t], (63) 

and define 9 and 9 S as, 

p (t) = P(0,t) ~r e °, t^+oo, (64) 
p s (t) = lim t0 _ +oo P(to, t) ~ r 0s , t -> +oo. (65) 

Po measures the first passage exponent of the growing interface, whereas p s contains the 
relevant information, when the interface has entered the steady state (to —> +oo). 

The correlators associated with these two persistence problems are respectively (when 
expressed as functions of the fictitious time r) [|17|j : 

f (T) = cosh(r/2f - | sinh(r/2)|^ (66) 
f s (T) = cosh^r/2) - ||2sinh(r/ 2)1", (67) 

and both satisfy 1 — /o,s(t) ~ for small r, with \i = 1 — d/z (// = 1 — (d + 2)/z for a 
volume conserving noise). We now simply treat /i as a free parameter. Using a connection 
to the fractional Brownian walker, it has been conjectured that 9 S = 1 — /i/2 ||17|| ) which has 
been confirmed by numerical simulations. 

Let us take two typical values for /i. For \x = 3/2, we find #o,var — 0.2088... and 
^o,pert — 0.2146..., which compare well to the numerical value 8q = 0.201 ± 0.005. For 
the case of the steady interface, the conjectured persistence exponent is 9 S = 1/4, in good 
agreement with the simulations (6 S = 0.247 ± 0.005). Variational and perturbative methods 
are reasonably accurate, giving S)Var = 0.2583 . . . and 9 SiPert = 0.2644 . . .. Note that the first 
order perturbation expression of Eq. (pg|) reproduces exactly the conjectured value for 9 S . 

We have also tested the case fj, — 3/4, which is getting dangerously close to the limit of the 
validity domain of our variational (yU var = 1/2) and perturbative (/i pe rt = 2/3) expressions. 
The numerical value of #o is 9q = 0.85±0.01, for which the variational approach gives #o,var = 
0.8852 . . .. Not surprizingly, the resummed perturbation leads to a bad result (#o,pert ~ 1-1). 
The conjectured value for 9 S is 9 S = 0.625, while the simulation of the process leads to 
9 S = 0.625 ± 0.005, and that of the associated discrete solid-on-solid model leads to 9 S = 
0.635 ±0.005 (see |T7| for details). We find a qualitatively correct value of S)Var = 0.6662 . . ., 
but the resummed perturbation fails again (0 SjPer t ~ 0.84). 

C. Smooth processes 

For singular processes (// < 2), it was not possible to compare our variational and 
perturbative results to the IIA expressions of Eq. (|4!|]44|), which are only defined for smooth 
processes. 
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One of the most spectacular example of smooth Gaussian processes has been given in 
12 1 : consider an initially random spatial distribution of charges of zero average, p(x, t = 0). 
It then evolves according to the simple diffusion equation, 

^(x,t) = V 2 p(x,t). (68) 

The persistence is defined as the probability that the local charge at a given x never changes 
sign. It decays as a power-law, defining 9d, the dimension-dependent persistence exponent. 
The IIA Jl2| (as well as the specific method of fL4| ) is in amazing agreement with numer- 



ical simulations. For instance, in d = 1, 9ua = 0.1203 . . ., to be compared to the numerical 
value 9 = 0.1207 ± 0.0005. The agreement seems to be of the same order in any dimension. 
Smooth processes are in principle beyond the range of application of perturbative methods. 
Still, the variational approach remains qualitatively correct, leading to 9 var = 0.1428. . . for 
the one-dimensional diffusion equation, whereas the resummed perturbation theory is again 
quite bad (# pe rt = 0.1612 . . .). 

Another example of smooth process is the Gaussian walker satisfying = r)(t), for 
n > 2 (n = 1 being the Markovian Brownian walker which is singular). The case n = 2 
corresponds to a particle submitted to a random force, for which the persistence exponent 
is known to be 9 = 1/4 p9| . The two-time correlator when expressed in the fictitious time 
reads, 

/(r) = ^exp(-|r|/2) - I exp(-3|r|/2). (69) 
The IIA leads to 9 nA = 0.2647. . ., whereas fl var = 0.2857. . ., and 9 pcrt = 0.3198 . . .. 



D. Perturbation around the action of a particle in a box 



Let us briefly give a few applications of our expressions of Eq. (|55| , ^7| , |59D . 

As a simple test, they have been applied to the case of the Markovian walker with 
A = 9 = 1/2, and Xo = 0. The "variational" approach, which consists in taking the size of 
the box (or fc ) such that the first order perturbation vanishes leads to 26* var = 1.0074. . .. 
For the correlator given by Eq. flBTD, we found # var = 1.4323..., in fair agreement with 
simulations and perturbative approaches around a Markovian process. 

Finally, we have tested Eq. ([59]) in the case of the Brownian walker (for which A = 1/2), 
for X 7^ 0. In this case, it is known that 9 satisfies T>2e(Xo) = |p7| , where V 2 e is a 
parabolic cylinder function. For X = 1/3, (comparable to (X 2 ) = /(0) = 1, for A = 1/2), 
we found 6* var = 0.7032 . . ., to be compared to the exact value 9 = 0.6440 . . .. Note that the 



perturbative expression of Eq. (45) leads to 9 = 0.6330 



IX. CONCLUSION 

In this paper, we have stressed the importance of studying persistence for Gaussian sta- 
tionary processes, as the calculation of 9 for many physical systems can be often mapped on 
the persistence problem for this kind of process. We have then extended the perturbative 
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approach around a Markovian process, introduced in jl(J. We have obtained resummed per 



turbative expressions (Eq. fl34|,37,51|)) and a new self-consistent perturbative (or variational) 
expression for the persistence exponent (Eq. (pEUP). It seems that this variational result is 
more effective in reproducing numerical results, sometimes with impressive accuracy. We 
have also shown that all these expressions take a similar form as the alternative result of the 
IIA, which only applies to smooth processes. We have also given perturbative expressions 
for the Xo-level persistence exponent (Eq. ( |2D|J^5| )). Finally, we have shown that this type of 
perturbative approach is even more general, as the starting process around which we decide 
to perturb can be any Markovian process associated with a (preferably solvable) quantum 
problem. We have illustrated this point by explicitly deriving a variational expression for the 
Xo-level persistence exponent, when the starting quantum system is chosen to be a particle 
in a bounded box (Eq. ( |55| , |57| , |59D ). 

Finally, we conclude by pointing out that our perturbative and variational techniques 
have been useful in a wide variety of problems. This includes the calculation of the survival 
probability of a mobile particle in a fluctuating field 



3C 



and the calculation of global 
persistence exponent in critical spin systems (to compute the order e 2 = (A — d) 2 perturbative 
correction [pi), and for directed percolation f31|. 
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APPENDIX A 



We want to compute Cj = |(6|a;|j)| 2 , where \j) is the jth eigenstate of the harmonic 
oscillator with an infinite barrier at the origin. 

(j\x) = V2(2j + l\x) = V* = g ( z ) e -#. (70) 

W+H2j + 1)!^ 



The extra factor \[2 is due to the fact that (j\x) is only defined on the interval [0; +oo], but 
should still be normalized. One then finds, 

Cj = 7r22i(2j + l)!^' (71) 

where Ij = J +o ° x 2 e~ x2 H 2 j+i(x) dx can be readily calculated, using the properties H n+1 {x) = 
2xH n (x) — 2nH n ~i(x), H' n (x) = 2nH n ^i(x) and fo°°e~ x H n (x) dx = ff n _i(0). This yields, 

which finally gives the result of the main text. 
The Cj's satisfy the recursion relation, 

= (2j - l) 2 
Cj (2j + 3)(2j + 2)' 1 } 

which shows that the generating function f(x) = Y^^ c j x ^ satisfies the hypergeometric 
differential equation x(l — x)f" + — \f = 0. The (unique) solution with /(0) = Co = - 
is given by 

m = -F(-i-\>%*)> ( 74 ) 



7T V 2' 2'2 



which yields the identities £g> = fi 1 ) = I and T^= Q jcj = /' (1) = f. 
Now defining {/(r) as in Eq. fl26|) we get, 



1 A 2 +oc A 2 

-(-d 2 + A 2 ) 2 C/(r) = yECi«P(-2jA|r|) = yS[exp(-2A|r|)], (75) 



(76) 



where Cj = (4j 2 — l) 2 c-,- satisfies (using Eq. (|7B|)), 

Cj+i _ J + I 



We recognize the recursion relation obeyed by the coefficients of the series expansion of the 
function S(x) = ^(1 — x)~ 3 / 2 , which leads to the final result of Eq. (|2~8D . 
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APPENDIX B 



In this appendix, we write the general equation for Eq and Ex, when the perturbation 
theory is applied to S = Sq+SS. Sq is assumed to be the action associated with the quantum 
Hamiltonian H , with eigenenergies E\ and eigenstates The associated Hamiltonian with 
an infinite wall at the origin is Hx, with eigenenergies ii and eigenstates When the 
potential is symmetric with respect to X = 0, one simply has (x\l) = y/2(x\2l + 1), and 



A. Equation for Eq 

The lower "energy" Eq is a functional of the inverse correlator g(ui) = 1/ f(u): 

1 +oo 

Eq = so + lim — - g (u)) ■ (0\X(-u n )X(u n )\0), (77) 

/3-+oo 2(3 n=0 

where, 

k(u n ) = ^ _1 K) = (0\X(-u n )X(u n )\0), (78) 

+°° Orl,m 2 

= J^2cosu n rE^\(0\x\l)\ 2 e-^-^dT = £ -*-L (79) 

with di = Si — Eq, mi = |(0|x|/)| (m = due to the symmetry of the potential). Note that 
k is nothing more than the Fourier transform of the two-time correlation function of the 
position of the considered quantum particle. 

Then, transforming the sum over Matsubara frequencies into an integral, one obtains: 



1 f+°° ( 
Ztt Jo \ 



^--l) dx. (80) 
/(*) / 



Due to the sum rule I]£=T 2d\m\ = 1, valid for any Hamiltonian, the integrand tends to 
as x — > +oo. 

In the text, we have considered two examples for the starting quantum correlator k: 

• Harmonic oscillator of frequency A: 

Ei = X(l + ~), and di = XI; mi = (2A)~ 1 / 2 5; ) i; k{uj) = (u 2 + A 2 ) -1 , which directly leads 



to Eq. (|T7|), obtained in section IV by expanding the exact result of Eq. (jig) . 



Particle in a box of width 2b = 7r/ko: 

Ei = |&o(Z + l) 2 , and di = ^^1(1 + 2). After an elementary calculation involving the 
eigenstates of a particle in a box, we get, 

„ / 1 4 (^+!) ii \ 256 +^ a.b, 

mi = ^-^^WTW' fc(w) = ^gwf^- (81) 

with, 
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2(i + l) 2 1 

a , = - zu 2 ±; - b ,- = -(2j + l)(2j + 3). (82) 

J 2j + 1 4 2j + 3 4 2 / v j v ) 



We then recover the expression of Eq. fl55"D . 

Note that the following sum rules have been used in the main text: 

256 +2? 256 +2? 1 256+2?% 5 vr 2 



B. Equation for E\ 

The "energy" E\ is also a functional of the inverse correlator g(u) = 1/ f(uj): 

^ +oo 

£ x = e + lim — £(SK) - ffoK)) • (6|X(-w n )*K)|6), (84) 

/3-+oo 2/3 n=Q 



where, k = % has been defined in Eq. ([78]). Let us now introduce, 

= (6|X(-cu n )XK)|6>, (85) 

+oo 9^7^5,2 

= /^cos^rE^ |(6|x|Z)| 2 e-^-^dr = £ ~^f\- (86) 

Z=0 «Z + 

^ is the two-time correlator of the position of the quantum particle in the presence of the 
wall. As before, d\ = E\ — Eq, rh\ = \(0\x\l)\. Note that, contrary to the calculation for Eq, 
the / = contribution in the sum above is non-zero so that, strictly speaking, this term 
should be written 2irm 2 ) 5(uj n ). This term has been written under this form in the main text. 
Finally, the general expression of E% reads, 

1 r +oc ( 1 l \ . 
E 1= i + — , K(x) dx. (87) 



2ttJo \f(x) k(x) 
We can again make this result more explicit for both considered quantum systems: 

• Harmonic oscillator of frequency A: 

ii — X(2l + |), d\ — 2X1, and rhi = y/cl, which leads to the expressions of Eq. (p^-|25|) . 

• Particle in a box of width 2b = tr/ho: 

In this case, K is the Fourier transform of the two-time correlator of the position of a 
particle in a box of size b = n/2k Q . It is then clear using Eq. fl8~T|) that, 

K(«>) = ^*(w/4) + ^6 (u) , (88) 

the extra Dirac peak coming from the fact that the operator X now has a finite average, 
as the particle belongs to the interval [0, b]. This immediately leads to the formula of 
Eq. (0), using fc(0) - 15 ~^ 
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When the constraint is X > Xo, with Xo = f]b (as in section VII), the quantum 
particle now lives in a box of size (1 — r])b, and the expression for K is changed 
accordingly, leading to, 

k{u) = £^j£% (**4^) + ^(1 + V) 2 S H • (89) 
This immediately leads to the result of Eq. (|59|) . 
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